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Abstract. In this paper, we describe new methods to compute the radical (resp. real 
radical) of an ideal, assuming it complex (resp. real) variety is finte. The aim is to combine 
approaches for solving a system of polynomial equations with dual methods which involve 
moment matrices and semi-definite programming. While the border basis algorithms of 
1171 are efficient and numerically stable for computing complex roots, algorithms based on 
moment matrices 1121 allow the incorporation of additional polynomials, e.g., to restrict the 
computation to real roots or to eliminate multiple solutions. The proposed algorithm can 
be used to compute a border basis of the input ideal and, as opposed to other approaches, 
it can also compute the quotient structure of the (real) radical ideal directly, i.e., without 
prior algebraic techniques such as Grobner bases. It thus combines the strength of existing 
algorithms and provides a unified treatment for the computation of border bases for the 
ideal, the radical ideal and the real radical ideal. 



1. Introduction 

Many problems in mathematics and science can be reduced to the task of solving zero- 
dimensional systems of polynomials. Existing methods for this task often compute all (real and 
complex) roots. However, often only real solutions are significant and one needs to sieve out all 
complex solutions afterwards in a separate step. 

Typical approaches in this vein are the efficient homotopy continuation methods in the spirit 
of |21| , [19j , recursive intersection techniques using rational univariate representation ^9] in the 
spirit of Kronecker's work [llj . Grobner basis approaches using eigenvector computations or 
rational univariate representation [S], [TB], [HI chap. 4]. In the latter methods, emphasis is 
put on exact input and computation. Using a different approach, Mourrain and Trcbuchct 
[17j have proposed an efficient numerical algorithm that uses border bases and the concept of 
rewriting family. In particular, in the course of this algorithm, a distinguishing and remarkable 
feature is a careful selection strategy for monomials serving as candidates for elements in a basis 
of the quotient space ]K[x]// (if / C K[x] is the ideal generated by the polynomials defining 
the equations). As a result, at each iteration of the procedure, the candidate basis for the 
quotient space ]K[x]// contains only a small number of monomials (those associated with a 
certain rewriting family). Another nice feature of this approach (and in contrast with Grobner 
base approaches) is its robustness with respect to perturbation of coefficients in the original 
system. 

On the other hand, Lasserre et al. |12] have proposed an alternative numerical method, 
real algebraic in nature, to directly compute all real zeros without computing any complex 
zero. This approach uses well established semi-definite programming techniques and numerical 
linear algebra. Remarkably, all information needed is contained in the so-called quasi-Hankel 
moment matrix with rows and columns indexed by the canonical monomial basis of K[x]c;. Its 
entries depend on the polynomials generating the ideal / and the underlying geometry when 
this matrix is required to be positive semi-definite with maximum rank. A drawback of this 
approach is the potentially large size of the positive semi-definite moment matrices to handle 
in the course of the algorithm. Indeed, when the total degree is increased from d to d + 1, the 
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new moment matrix to consider has its rows and columns indexed by the canonical (monomial) 
basis of K[x]c(_|-i. 

The goal of this paper is to combine a main feature of the border basis algorithm of |17) 
(namely its careful selection of monomials, considered as candidates in a basis of the quotient 
space K[x]//) with the semi-definite approach of [12] for computing real zeros and an approach 
for computing the radical ideal inspired by jlOj . 

The main contribution of this paper is to describe a new algorithm which incorporates in the 
border basis algorithm the positive scmi-definitencss constraint of the moment matrix, which 
arc much easier to handle than the relaxation method of [12] . We show the termination of the 
computation in the case where the real radical is zero-dimensional (even in cases where the 
ideal is not zero-dimensional). A variant of the approach is also proposed, which yields a new 
algorithm to compute the (complex) radical for zero-dimensional ideals. 

In this new algorithm, the rows and columns involved in the semi-definite programming 
problem are associated with the family of monomials (candidates for being in a basis of the 
quotient space) and its border, i.e., a subset of monomials much smaller than the canonical 
(monomial) basis of E[x]ji considered in [T^. As a result, the (crucial) positive scmi-definiteness 
constraint is much easier to handle and solving problem instances of size much larger than those 
in [12| can now be envisioned. A preliminary implementation of this new algorithm validate 
experimentally these improvements on few benchmarks problems. 

The approach differs from previous techniques such as [I] which involve complex radical com- 
putation and factorisation or reduction to univariate polynomials, in that the new polynomials 
needed to describe the real radical arc computed directly from the input polynomials, using 
SDP techniques. 

The paper is organized as follows. Section [2] recalls the ingredients and properties involved 
in the algebraic computation. Section [3] describes duality tools and Hankel operators involved 
in the computation of (real) radical of ideals. In Section HJ we analyse the properties of the 
truncated Hankel operators. In section [S] we describe the real radical and radical algorithms 
and prove their correctness in section [SJ Finally, Section [7] contains some illustrative examples 
and experimentation results of a preliminary implementation. 

2. Polynomials, dual space and quotient algebra 

In this section, we set our notation and recall the eigenvalue techniques for solving polynomial 
equations and the border basis method. These results will be used for showing the termination 
of the radical border basis algorithm. 

2.1. Ideals and varieties. Let K[x] be the set of the polynomials in the variables x = 
(xi, . . . , x„), with coefficients in the field K. Hereafter, we will chooscQ K = R or C. Let 
IK denotes the algebraic closure of K. For a G N", x" = x"^ • • • a;"" is the monomial with expo- 
nent a and degree |a| = X)i '^i- The set of all monomials in x is denoted M. = (x). We say 
that x" < -xP if x" divides x^, i.e., if a < /3 coordinate- wise. For a polynomial / = /qX", its 
support is supp(/) := {x" | /„ ^ 0}, the set of monomials occurring with a nonzero coefficient 
in /. 

For < e N and S C K[x], we introduce the following sets: 

• St is the set of elements of S of degree < i, 

• S'[t] is the set of element of S of degree exactly t, 

• (S") = { X]/(E5 ^/ / I / S 5", A/ e K} is the linear span of 5, 

• (5) = { J2fesPf f \ Pf ^ IK[x], / e 5*} is the ideal in K[x] generated by S, 



For notational simplicity, wc will consider only these two fields in this paper, but R and C can be replaced 
respectively by any real closed field and any field containing its algebraic closure) 
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• {S\t) = {J2feStPf -f I Pf ^Wt-dcg(/)} is the vector space spanned by {x"/ | / S 
^t,|a| <t-deg(/)}, 

• 5"+ := S* U xiS U . . . U XnS is the prolongation of S by one degree, 

• dS ~S+\S is the border of 5, 

• := 5*^ is the rcsuh of applying t times the prolongation operator on S, with 
S-Iil = 5"+ and, by convention, = S. 

Therefore, St = SnK[x]t, S^] = Sr\K[x]it], S*'*' = {x^f | / e 5, |a| < t}, {S\t)C {S)r\K[x]t = 
{S)t, but the inclusion may be strict. 

li B C contains 1 then, for any monomial m E A4, there exists an integer k for which 
m € B^'^y The B-index of m, denoted by St3{m), is defined as the smallest integer k for which 

A set of monomials B is said to be connected 1 if 1 £ S and for every monomial m 1 in 
B, m ~ Xi^m! for some io G [1, n] and m! € B. 

Given a vector space E C K[x], its prolongation := E + xiE + . . . + XnE is again a vector 
space. 

The vector space E is said to be connected to 1 if 1 G -B and any non-constant polynomial 
p G E can be written a.s p = po + J27=i ^iPi ^'^^ some polynomials poiPi S E with deg(po) < 
deg(p), degijji) < deg(p) — 1 for i € [1, n]. Obviously, £^ is connected to 1 when E = (C) for 
some monomial set C C which is connected to 1. Moreover, i?+ = (C+) if i? = (C). 

Given an ideal / C K[x] and a field L 3 K, we denote by 

:= {z e L" I /(x-) - V/ e /} 

its associated variety in L". By convention V{I) = For a set F C K", we define its 

vanishing ideal 

I{V) := {/ e K[x] I f{v) = Vi; e F}. 

Furthermore, we denote by 

VI := {/ e K[x] I /™ G / for some m £ N \ {0}} 

the radical of /. 

For K = M, we have V{I) = Vc{I), but one may also be interested in the subset of real 
solutions, namely the real variety Vr(/) = y{I) H M". The corresponding vanishing ideal is 
^(^(^)) and the real radical ideal is 

\/l — {pe M[x] I + ^ g2 g J fQj, gQj^g g m e N \ {0}}. 

Obviously, 

/c VJc /(yc(/)), /c /(14(/)). 

An ideal / is said to be radical (resp.. real radical) li I — \fl (resp. / = v^). Obviously, 
/ C /(!/(/)) C Hence, if / C R[x] is real radical, then / is radical and moreover, 

V{I) = VvXl) C M" if < oo. 

The following two famous theorems relate vanishing and radical ideals: 

Theorem 2.1. 

(i) Hilbert's Nullstellensatz (see, e.g., [6j §4.1]; \fl = /(VcC^)) /or an ideal I C C[x]. 

(ii) Real Nullstellensatz (see, e.g., [3l §4.1]^ \fl = /(Fr(/)) /or an ideal I C R[x]. 
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2.2. The quotient algebra. Given an ideal / C ]K[x], the quotient set IK[x]// consists of all 
cosets [f] ■= f + 1 = {f + q \ q G /} for / G IK[x], i.e., all equivalent classes of polynomials of 
K[x] modulo the ideal /. The quotient set IK[x]// is an algebra with addition [/] + [g] := [f + g], 
scalar multiplication A[/] :~ [A/] and with multiplication [f][g] := [fg], for A G M, /, g G K[x]. 

A useful property is that, when / is zero-dimensional (i.e., < oo) then K[x]// is 

a finite-dimensional vector space and its dimension is related to the cardinality of V{I), as 
indicated in Theorem 12.21 below. 

Theorem 2.2. Let I be an ideal in K[x]. Then |Vj^(/)| < oo dimK[x]// < oo. Moreover, 
|^(-^)| l£ dim ]K[x]//, with equality if and only if I is radical. 

A proof of this theorem and a detailed treatment of the quotient algebra K[x]// can be found 
e.g., in [6], E], m- 

Assume |Vj^(/)| < oo and set N := dimK[x]//, |Vj^(/)| < TV < oo. Consider a set B := 
{bi, . . . , 6jv} ^ IK[x] for which {[bi], . . . , [few]} is a basis of K[x]//; by abuse of language we also 
say that B itself is a basis of K[x]//. Then every / G K[x] can be written in a unique way 
as / = '^i^i + ?! where c; G K, p G /; the polynomial ni^B^f) J2i=i '^i^i called the 

remainder of / modulo /, or its normal form, with respect to the basis B. In other words, (B) 
and K[x]// are isomorphic vector spaces. 

2.2.1. Multiplication operators. Given a polynomial h G K[x], we can define the multiplication 
(by h ) operator as 

Xh-. IK[x]// K[x]// 
^ ' [/] ^ X,{[f]) [hf], 

Assume that N :— dimK[x]// < oo. Then the multiplication operator can be represented 
by its matrix, again denoted Xh for simplicity, with respect to a given basis B = {&i, . . . , b^} 
of K[x]//. 

Namely, setting 7r/_g(ft,5j) ^iLi'^ij^i for some scalars Oij G K, the jth column of X^ 
is the vector {aij)fLi. Define the vector (is^y ibj{v))jLi G K^, whose coordinates are the 
evaluations of the polynomials bj G B at the point w G K . The following famous result (see 
e.g., O Chapter 2§4], [8]) relates the eigenvalues of the multiplication operators in IK[x]// to 
the algebraic variety V{I). This result underlies the so-called eigenvalue method for solving 
polynomial equations and plays a central role in many algorithms, also in the present paper. 

Theorem 2.3. Let I be a zero-dimensional ideal in ]K[x], B a basis o/]K[x]//, and h G K[x]. 
The eigenvalues of the multiplication operator Xh are the evaluations h{v) of the polynomial h 
at the points v G V{I). Moreover, {Xh)'^CB,v ~ h[v)(^B.v and the set of common eigenvectors of 
{Xh)heK[x] are up to a non-zero scalar multiple the vectors (^g.^, for v G V{L). 

Throughout the paper we also denote by Xi := Xx^ the matrix of the multiplication operator 
by the variable Xi. By the above theorem, the eigenvalues of the matrices Xi are the ith 
coordinates of the points v G V{I). Thus the task of solving a system of polynomial equations 
is reduced to a task of numerical linear algebra once a basis of K[x]// and a normal form 
algorithm are available, permitting the construction of the multiplication matrices Xi. 

2.3. Border bases. The eigenvalue method for solving polynomial equations from the above 
section requires the knowledge of a basis of K[x]// and an algorithm to compute the normal 
form of a polynomial with respect to this basis. In this section we will recall a general method 
for obtaining such a basis and a method to reduce polynomials to their normal form. 

Throughout S C is a finite set of monomials. 

Definition 2.4. A rewriting family F for a (monomial) set B is a set of polynomials F = 
{fi}iei such that 
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. supp(/,) C B+, 

• /, has exactly one monomial in dB, denoted as j{fi) and called the leading monomial 
of fi- (The polynomial fi is normalized so that the coefficient of^{fi) is 1.) 

• if lift) = 7(/j) then i = j. 

Definition 2.5. We say that the rewriting family F is graded if dcg{-f{f)) = deg(/) for all 
f^F. 

Definition 2.6. A rewriting family F for B is said to be complete in degree t if it is graded 
and satisfies (dB)t C ^{F); that is, each monomial m S dB of degree at most t is the leading 
monomial of some (necessarily unique) f €z F. 

Definition 2.7. Let F be a rewriting family for B, complete in degree t. Let t:f,I3 be the 
projection on (B) along F defined recursively on the monomials m ^ Ait in the following way: 

• if m €z Bt, then 7ri?.B(m) = m, 

• ifm e {dB)t {B^^\B^^^)t), then ttf,b{iti) = where f is the (unique) polynomial 
in F for which j{f) ~ m, 

• if me (fiW \6['=-il)t for some integer k > 2, write m = XigTn' , where m' £ ' and 
io e is the smallest possible variable index for which such a decomposition exists, 
then 7rF,e(m) = nF.Bixi,, nF,Bini')). 

One can easily verify that deg{TTF^B{^)) ^ deg(m) for m G Mt- The map ttf,b extends by 
linearity to a linear map from K[x]t onto {B)f By construction, / = j{f) — iTFfiilif)) and 
'^F,B{f) = for all / g Ft- The next theorems show that, under some natural commutativity 
condition, the map 'nF,B coincides with the linear projection from K[x]t onto {B)t along the 
vector space {F\t), and they introduce the notion of border bases. 

Definition 2.8. Let B C A4 be connected to I. A family F C K[x] is a border basis for B if it 
is a rewriting family for B, complete in all degrees, and such that K[x] = (B) ® (F). 

An algorithmic way to check that we have a border basis is based on the following result, 
that we recall from [17] : 

Theorem 2.9. Assume that B is connected to 1 and let F be a rewriting family for B , complete 
in degree t € N. Suppose that, for all m £ Mt-2, 

(2) 7ri?3(xi 7ri?.g(2;j to)) = 7ri?,g(a;j 7ri?.g(xi m)) for all i, j E [l,n]. 

Then ttf.b coincides with the linear projection o/]K[x]i on {B)t along the vector space {F\t); 
that is, IK[x]t = {B)t® {F\t). 

Proof. Equation ([2]) implies that any choice of decomposition of m E Ait as a product of 
variables yields the same result after applying ttf.b- Indeed, let to ~ Xi-^ m' = Xi^ m" with ii ^ 
12 and m',m" £ Ait-i- Then there exists to'" e Ait-2 such that m' = Xi^ to'", m" = Xi^ m!" . 
By the relation ([2]) we have: 

T^FfiiXii TTF,B[m')) 

= T^Ffiixii T^Ffiixi^m'")) = ttf,b{xi^ T^Ffiixi^m'")) 
= t^f,b{xi^ T^Ffiim")). 
Let us prove by induction on I = deg(m) that for a monomial m = Xi^ ■ ■ ■ Xi^ € A4t, 

(3) 7ri^,e(™) = TTF,BixnTTF,Bixn ■ ■ ■ T^FMix^) ■■ ■ ), 

does not depend on the order in which we take the monomials in the decomposition m = 

•^ii ' ' ' '^ii • 

• Either to E B. As S is connected to 1, there exists i' E [1, n] and to' E St-i such that 
TO = ttf,b(j^) = T^F,B{xi'fTi') — T^F.Bi^i' ttf.sIti')), from which we deduce (j3]) using the 
induction hypothesis applied to m' and relation 
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• Or m ^ B. Then, by definition of ttp^b, there exists i' e [1,?^] and m' G Mt~i such 
that 7ri?.g(7Ti) = 7ri?.e(-'^i' i'F.e(™'))i from which we deduce ([3]) in a similar way using 
the induction hypothesis appHed to m' and relation ([2|). 

The map ttf,b defines a projection of K[x]t on {B)t. It suffices now to show that Kcr7ri?.g = 
{F\t). First we show that m—TTp^js{rn) S {F\s) for all m G A^^, using induction on s = 0, . . . , i. 
The base case s = is obvious; indeed TrF,B{'^) = 1 since 1 E B, and G (F|0). Consider 
TO € Ms+i- Write to = Xi^m! where to' € and ttf,b{'>^'>') ~ '^F.eixio'^F.Bi'm-')) (recall 
Definition 12. 7p . We have: 

TO - TTF,B{m) = Xi„{m' - TTF.Bim')) + XigTlF,B{m') - T^F.B{XiaTTF.B{m')) . 

^ ^ ^ ^ ^ ^ 

■.=q ■.=!• 

By the induction assumption, m' —7Tf.b{itt-') G | s) and thus q G {F \ s + 1). Write TTF.Bini') = 
Shes, "^bb [Xb e K). Then, r = J2beBs ^bixigb - TTF,B{xiob)), where Xi^b - TTF,B{xi„b) = if 
Xigb € B, and Xi^^b — T^F.B{xi„b) is a polynomial of i^s+i otherwise. This implies r G {F \ s + \) 
and thus to — T^F,B{Tn) G (i^ | s + 1). Thus we have shown that IK[x]t = {B)t + {F\t). Next, 
observe that {F\t) C KerTTi^^g, which follows from the fact that Ft C Y^ei-KF.B together with 
©. This implies that {B)t n{F\t) = {0} and thus the equality {F\t) = Kqi-kf^- □ 

In order to have a simple test and effective way to test the commutation relations we 
introduce now the commutation polynomials. 

Definition 2.10. Let F he a rewriting family and f,f' G F. Let m,m' be the smallest de- 
gree monomials for which mj{f) ~ m'j{f'). Then the polynomial C{f,f') := mf — m' f = 
m''KF,B{f') ~ 'm-'^F,Bif) called the commutation polynomial of f,f'. 

Definition 2.11. For a rewriting family F with respecet to B, we denote by C^{F) the set of 
polynomials of the form m f — m' /', where f, J' G F and to, to' G {0, 1, Xi, . . . , x„} satisfy 

• either m^{f) ~ m' "f{f'), 

• or m"/{f) G B and to' = 0. 

Therefore, C+(F) C {B^) and C^{F) contains all commutation polynomials C{f,f') for 
f,f'€F whose monomial multipliers to, to' are of degree < 1. The next result can be deduced 
using Theorem 12.91 

Theorem 2.12. Let B d Ai be connected to 1 and let F be a rewriting family for B, complete 
in degree t. If for all c G C^{F) of degree < t, TrF,B{c) = 0, then ttf.b is the projection o/]K[x]t 
on {B)t along {F\t), ie. ]K[x]t = {B)t ® (F 1 1). ' 

Proof. Let us prove by induction on t that if F is complete in degree t and for all c G C+ [F] of 
degree < t, 'KF.Bic) = then any to G M.t-2 satisfies which in view of Theorem 12.91 suffices 
to prove the theorem. 

Let us first prove that ^ holds for to G Bt-2- We distinguish several cases. If x^to, Xjm G B 
then ^ holds trivially. Suppose next that Xim,Xjm G dB. Then, / :— Xim — i^F.Bixim) and 
/' := Xjm — -Kpfiixjm) belong to Ft-i- As Xj^{f) — Xij{f'), xjf — Xif G C~^{F) and thus, 
by our assumption, ■nF,B{xjf) = T^F.B{xif'), which gives ([5]). Suppose now that x^to G dB and 
Xjm G B. As before / = .t,;to — TTF,Bixim) G Ft-i. If Xjj{f) = XiXjm G B then Xjf G C^{F) 
and thus TTF,B{xjf) = gives ([5]). Otherwise, xtXjm G dB and let /' := XiXjm~TrF,BixiXjm) G 
Ft. Now, Xjj{f) = 7(/') implies Xjf — /' G C~^{F) and thus TTF,B{xjf — /') = which gives 
again ([2]). This shows ([2|) in the case when to G Bt-2, and thus we have 

(4) '^F.Bixi2'^F.Bixiib)) ^ TTFfiixi^T^F.Bixi^b)) for all b G {B)t-2- 

Let us now consider to G A4t~2\Bt~2- By definition T:F,B{xi rn) ~ T^F^sixi' T^F^BifTT-')) for 
some to' G and i' G [1, n] such that Xim = Xi' to'. If i 7^ i' there exists to" G such 
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that m = Xi' m", m' = Xi to". As F is also complete in degree t — 1 and for all c G C+(F) of 
degree < t — 1, 'n'F,B{c) = 0, by induction hypothesis we have 

so that 7rir_B(a;i to) = HF.Bi^i'^F.BifTi)). If i = «', we have by definition 7ri?.g(a;i to) = '^F.Bixi'^F,Bi' 
As ■KF.B{fTi) = fTi for TO G Bt~2, we deduce that 

(5) '^F,Bixim) ~ iTF,B{xi''^F,B{'m)) for all rii £ i G [I,"-]- 

Now, using ([5]), 'KF.Bi^iT^F.Bi^j'm)) is equal to 7ri?_B(a;i7ri?_B(a:;j7ri? ^(to))) which in turn is 
equal to tt^? ^(xjTrp g(xj7ri? ^(to))) (using (|3])) and thus to F^BixjT^ F.Bixim)) (using again ([5])). 
We can now apply Theorem 1 2 . 91 and conclude the proof. □ 

Theorem 2.13 (border basis, [T7]). Let B C A4 be connected to 1 and let F he a rewriting 
family for B, complete in any degree. Assume that TTF,Bic) = for all c € C'^{F). Then B is 
a basis o/K[x]/(F), K[x] = (B) ® (F), and {F)t = {F\t) for all t G N; the set F is a border 
basis of the ideal I = [F) with respect to B. 

Proof. By TheoremHH ]K[x]t = {B)t®{F\ t) for aU t e N. This imphcs that K[x] = (B) ® (F) 
and thus S is a basis of K[x]/(F). Let us prove that {F)t = {F\t) for all t £ N. Obviously, 
{F\t) C {F)t. Conversely let p <E {F)t. Then p = r + q, where r G {B)t and q & {F \ t). Thus 
p - g e (F) n (6) = {0}, i.e., p^q£{F\t). □ 

This implies the following characterization of border bases using the commutation property. 

Corollary 2.14 (border basis, |16]). Let B d Ai be connected to 1 and let F be a rewriting 
family for B, complete in any degree. If for all m £ B and all indices i,j G we have: 

TTF,B{Xi T^F,B{Xj to)) = TTf,b{Xj TTf,b[Xi to)), 

then B is a basis o/K[x]/(F), K[x] = {B) ® [F), and {F)t = {F\t) for all t e N. 

Proof. Same proof as for Theorem 12. 13i using Theorem 12.91 □ 

3. Hankel Operators 

In this section, we analyse the properties of Hankel operators and related moment matrices, 
that wc will need hereafter, for the moment matrix approach. 

3.1. Linear forms on the polynomial ring. The set of K- linear forms from K[x] to K is 
denoted by K[x]* := HomK(K[x],K) and called the dual space of K[x]. A typical element of 
K[x]* is the evaluation at a point ( £ K": 

1^ : pe K[x] ^ p{C) e K. 

Such evaluation can be composed with differentiation. Namely, for a E N", the differential 
functional: 

evaluates at ( the derivative d" of p. For a = 0, ■ d'^ = The dual basis of the monomial 
basis (x")QgN" of K[x] is denoted (d")Q,gNn; we have d"(x^) = Sa.p- In characteristic 0, 
d" := lo • j-j„^ a ! ^"- Any element A e IK[x]* can be written as A = J2aM^°')<i°' ■ In 
particular, 1^ = EaeN" C"d". 
For S' C K[x], we define 

:= {A e IK[x]* \ ypeS A{p) = 0}. 
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3.2. Hankel operators. The dual space IfC[x]* has a natural structure of K[x]-niodule which 
is defined as follows: (p, A) e K[x] x IK[x]* p • A e IK[x]*, where 

p • A : q e K[x] ^ A{pq) G K. 

Note that, for any a, (3 € N", we have 

x^-d" = d"-'' ifa>^, 

= otherwise . 

Definition 3.1. For A G ]K[x]*, the Hankel operator Ha is the operator from IK[x] to K[x]* 
defined by 

Ha : peK[x]^ p- AeK[x]*. 

Lemma 3.2. For A G K[x]*, the matrix of the Hankel operator Ha with respect to the bases 
(x") o/K[x] and (d'^) ofK[x]* is [Ha] = (A(x"+^)). 

Proof. Writing A = J^-y A(x''')d'''. we have: 

i7A(x")=x"-A = ^A(x'')x"-d''= J2 A(x'^)d''-" =^A(x"+^)d''. 

7 7l7>a 13 

□ 

We now summarize some well known properties of the kernel 

KevHA = {pe ]K[x] I p • A = 0, i.e., A{pq) = Vg G K[x]}. 

of the Hankel operator Ha- Recall the definition of a Gorcnstcin algebra [3], [51 Chap. 8]. 

Definition 3.3. An algebra A is called Gorcnstcin if A and its dual space A* are isomorphic 
A-modules. 

Applying this definition to ^ := ]K[x]/Kcri/A yields 

Lemma 3.4. Kcr_ffA is an ideal in K[x] and the quotient space A := K[x]/Ker_ffA o Gorcn- 
stcin algebra. 

Proof. Direct verification, using Ha as isomorphism in the proof of the second part of the 

lemma. □ 

The focus of this paper is the computation of zero-dimensional varieties, which relates to 
finite rank Hankel operators as shown in the following lemma. 

Lemma 3.5. The rank of the operator Ha is finite if and only i/Ker TJa is a zero- dimensional 
ideal, in which case dimK[x]/ Ker_ffA = rankiJA- 

Proof. Directly from the fact that, given pi, . . . ,Pr € IK[x], Ha{pi), • • • , -ffA(Pr) are linearly inde- 
pendent in K[x]* if and only if the cosets [pi], . . . , [pr] are linearly independent in IK[x]/ Kcr/^A- 

□ 

The next theorem states a fundamental result in commutative algebra, namely that all zero- 
dimensional polynomial ideals can be characterized using differential operators (see [51 Chap. 7], 
[H Thm. 2.2.7]). For the special case of zero-dimensional Gorcnstcin ideals, a single differential 
form is enough to characterize the ideal. 

Tiieorem 3.6. Let K = C and assume rank Ha = r < oo. Then there exist d, . . . G C" 
(with d < r) and non-zero (differential) polynomials pi,...,pd G C[9], of the form Pi{d) = 
X^qgA ^i,ad°' where Ai C N" is finite and a.i^a G such that 

d 

(6) A = 5]1c.-k(9). 

i=l 
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For a zero-dimensional ideal / C ]K[x] with simple zeros V{I) = {Ci, ■ • ■ , Cr} C K" only, we 
have = {Iq-^ , . . . , 1^^) and the ideal / is radical as a consequence of Hilbert's Nullstellensatz. 

In a similar way, we can now characterize the linear forms A for which Keri?A is a radical 
ideal. 

Proposition 3.7. Let K = C and assume that ranki^A = r < oo. Then, the ideal Ker H\ is 
radical if and only if 

r 

(7) ^ = Y1 ^^^^ Aj e K - {0} and Q G IK" pairwise distinct, 

1=1 

in which case KeriJA = ^(Ci? • ■ • i Cr) the vanishing ideal of the Ci s- 

Proof. Assume first that Keri^A is radical with V^(Keri7A) {Ci: • • ■ : Cr} C IK". This implies 
KeriJA = I{V{Kcr H\)) = I{Ci,...,(r)- Let pi £ C[x] be interpolation polynomials at the 
points (i, i.e., Pi{Cj) = foi' h.i ^ Then the set {pi, . . . ,pr} is linearly independent in 
A := K[x]/(KeriJA) and thus is a basis of A. As the linear functionals A and X)I=i ^(pO^Ci 
take the same values at each pi, we obtain: A = ^{Pi)'^Ci- Moreover, := A{pi) ^ 0, 

since ranki^A = ^■ 

Conversely assume that A is as in ([7]). The inclusion I{Ci,...,Cr) C Ker_ffA is obvious. 
Consider now p g Ker H\ and as before let Pi G IK[x] be interpolation polynomials at the 
(i's. Then = A{ppi) = Kp{Ci) implies p{(i) = 0, thus showing p £ /(Ci, . . . , Cr)- As 
KeriJA = -^(Cii ■ • • 7 Cr) is the vanishing ideal of a set of r points, it is radical by the Hilbert 
Nullstellensatz. □ 

In a similar way, we can also characterize real radical ideals using Hankel operators. 

Proposition 3.8. Let K = M and assume that rankiJA — r < oo. Then, the ideal KeriJA is 
real radical if and only if 

r 

(8) A = ^Aac, with X, e R - {0} and e pairwise distinct. 

i=l 

Proof If Keri?A is real radical then V^(Keri7A) = {Ci, . • . , Cr} C M", so that ^ gives ©. Con- 
versely, if A is as in dS]), then Kei-H\ is real radical, since J2j l] ^ Ker if a implies J2j 9j(Cj)^ = 
and thus qjiCi) = 0, giving qj G KcriJA- D 

Let us now recall a direct way to compute the radical of the ideal Ker H\. First, consider 
the quadratic form Qa defined on K[x] by 

(9) Qa : {p,q)eK[^]^^A{pq)eK. 

Then, QA{p,q) — A{pq) ~ HA{p){q) = HA{q){p) for all p,q & IK[x], and the matrix of Qa 
in the monomial basis (x") is [Qa] = (A(x"+'^)). We saw in Lemma 13.41 that the algebra 
A = K[x]/Ker7?A is Gorenstein. An alternative characterisation of Gorenstein algebras states 
that the above quadratic form Qa defines a non-degenerate inner product on A (see eg. [5] [chap. 
9]). Assume now that rankiJA = r < oo so that dim„4 = r. Let &i, . . . , &r be a basis of A and 
let di, . . . , d,- be its dual basis in A for Qa: it satisfies A(6i dj) = 5ij for i,j G [1, r]. Then, for 
any element a A, we have 

r 

(10) a = ^A{ad,)b,. 

i=l 

In particular, we have the following property: 

Proposition 3.9. Let A := X]I=i ^i - Given h G A, let Xh he the corresponding multiplication 
operator in A. We have 

Trace(A'ft) = A(/iA). 
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Proof. By relation ([T(I|. the matrix of Xh in the basis {bi)i<i<r of A is {A{hbj di))i<ij<r and 
thus its trace is 

r 

Trace(A',,) = ^ A(/i&,dO = A{hA). 

i=l 

□ 

As a direct consequence we deduce the following result (see e.g., [TU]): 

Theorem 3.10. Let IK = C and assume rank7?A ~ r < oo. Let fei, . . . ,6^ be a basis of Aa, 
di, . . . ,dr be its dual basis with respect to the inner product given by Q\, and A = X]i=i ^i- 
Then the radical of¥^eTH\ is Kev H^.\. 

Proof. Let / := KctHa. A polynomial h is in y/l if and only if some power of h is in I 
or, equivalently, if and only if Xh is nilpotent. By a classical algebraic property, the latter is 
equivalent to Trace(A';iA'a) = = Tracc(A'/ia) for all a £ A. Indeed, as the operators Xh,Xa 
commute, if Xh is nilpotent then so is X^Xa and we have TT:&ce{XhXa) = 0. Conversely if 
Ti:ace{XhXa) = for all a G „4 then, by Cayley-Hamilton identity, the characteristic polynomial 
det(A/ — Xh) of Xh is A*" and thus Xh is nilpotent. By Proposition [2111 we deduce that h £ \fl 
if and only if A(A ha) ~ for all a e ^a, that is, if and only if e Ker_ffA A- D 

3.3. Positive linear forms. We now assume that IK = R and consider the polynomial ring 
M[x]. We first show that the kernel of a Hankel operator Ha is a real radical ideal when 
A g R[x]* is positive. This result is crucial in the algorithm that computes the real radical of 
an ideal. 

Definition 3.11. We say that A G M[x]* is positive, which we denote A > 0, if A{p'^) ^ for 
all p G R[x]. Equivalently, we will say Ha ^ i/ A ^ 0. 

We will use the following simple observation. 

Lemma 3.12. Assume A )>= 0. For p G R[x], A(p^) = implies p G KcriJA and thus A{p) = 0. 
For A, A' :>= 0, Kcr Ha+A' = Kcr Ha n Ker Ha' ■ 

Proof For any q G R[x], t G R, A{{p + tqf) = t^ A{q^) + 2tA{p q) > 0. Dividing by t and letting 
t go to zero yields A{pq) = 0, thus showing p G KeriJA- The inclusion Ker Ha H KbiHa' C 
KerTJA+A' is obvious. Conversely, let p G KerTJA+A'- In particular, (A + A')(p^) = 0, which 
implies A{p'^) = A'(p2) = (since A{p'^) , A' {p'^) > 0) and thus p G KcriJA n Kcri^A'- □ 

Proposition 3.13. // A 0, then Kcr Ha is o, real radical ideal. 

Proof. Assume J^iPf ^ Kcr Ha; we show that pi G Ker Ha. Indeed, {J^iPi) ■ A = implies, 

for all q G R[x], = A{Y,,piq^) = T,i HpIq^) and thus A{pjq^) = 0. By LemmaEia this in 

turn implies A{piq) = and thus pi G Ker Ha- O 

We saw in Proposition |3lH] that the kernel of a finite rank Hankel operator Ha is real radical 
if and only if A is a linear combination of evaluations at real points. We next observe that A is 
positive precisely when A is a conic combination of evaluations at real points. 

Proposition 3.14. Assume rank Ha ~ r < oo. Then A)p Q if and only if A has a decomposi- 
tion (0) with \i > and distinct Q G R", in which case V (Ker Ha) = {Ci, • ■ • , Cr} C R". 

Proof. If A = Yli=i ^i'^Ci with > and Q G R", then A )== holds obviously. Conversely, 
assume that A ^ then by Proposition l3.13l the ideal Ker Ha is real radical. By Proposition l3.8| 
A has a decomposition ^ where Xi ^ A{pi) 0, Q £ R", and pi arc interpolation polynomials 
at the Ci's. As p^ — pi G /(Ci, ■ • ■ , Cr) = Keri^A, we have A{pi) = A(p?) > 0, which concludes 
the proof. □ 
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To motivate the next section, let us recall Lemma 13.51 and observe how it specializes to 
truncated Hankel operators defined on subspaces of K[x]: 

Lemma 3.15. Let B ~ {bi, . . . , br} C M[x] and A G ]K[x]*. The operator 

: {B) ^ {BY 

r 

P = ^ Xibi ^ p - K 

i=l 

has a trivial kernel if and only if the cosets . . . , [br] G IK[x]/ KeriJA are linearly independent 
m IK[x]/KerHA. 

Proof. Direct verification using the fact that Kcr = Kcr Ha n (6) . □ 

Assume now that Kcr Ha is zero-dimensional and that B = {bi,.. .,br} C K[x] is chosen 
so that [bi], . . . , [br] form a basis of ^ = K[x]/Keri7A- As in relation (|9]), we consider the 
quadratic form on A defined by 

Qa ■ {p,q) e Ax A{pq) e K. 

Note that a matrix representation of this form can be obtained by taking the principle submatrix 
of [Ha] indexed by B. Following [Tl], we recall under which conditions the bilinear form 
relates to the Hermite form 

Th-.Ax A^K 

if,g) ^ TraceiXf gh) 

for some h E A. 

Lemma 3.16. The quadratic form associated to coincides with the Hermite form Th for 
some h € A if and only if Kcr Ha is radical. 

Proof See [H Sec. 2.2]. □ 

4. Truncated Hankel Operators 

We have seen in the previous section that the kernel of the Hankel operator associated to a 
positive linear form is a real radical ideal. However, in order to be able to exploit this property 
into an algorithm, we need to restrict our analysis to matrices of finite size. For this reason, we 
consider here truncated Hankel operators, which will play a central role for the construction of 
(real) radical ideals. 

For E C K[x], set E ■ E {pq [ p,q €z E}. Suppose now E C K[x] is a vector space. A 
linear form A defined on {E ■ E) yields the map H^ : E ^ E* hy Hf (p) = p • A for p G E. 
Thus Hf^ can be seen as a truncated Hankel operator, defined only on the subspace E. 

Given a subspace Eq C E, A induces a linear map on {Eq ■ Eq) and we can consider the 
induced truncated Hankel operator Hj^" : Eq — > (£^o)*- 

Definition 4.1. Given vector subspaces Eq C E C K[x] and A G {E ■ E)* , H^ is said to be a 
flat extension of its restriction H^" to Eq if rank H^ = rank H^° . 

We now give some conditions ensuring that it is possible to construct a flat extension of a 
given truncated Hankel operator. The next result extends an earlier result of Curto-Fialkow 
[7|; a generalization of this result can be found in [5]. 

Theorem 4.2. [15j Consider a vector subspace E C K[x] and a linear function A on {E^ ■ E^) . 

Assume that E = (C) where C CZ M is connected to 1 and that r&nkHf = rankiJ^. Then 
there exists a (unique) linear function A G ]K[x]* which extends A, i.e., A{p) = A(p) for all 
p G {E^ '^'^'^ satisfying rank Hj^ = rankH^^ . In other words, the truncated Hankel 

operator H^^ has a (unique) flat extension to a (full) Hankel operator Hj^. 
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In the following, we will deal with linear forms vanishing on a given set G of polynomials. 
Definition 4.3. Given a vector space E C ]K[x] and G d {E ■ E), define the set 

(11) Cg,e {Ae{E- EY \ A(g) = e G}. 
// IK = M, define 

(12) Cg,e,> ■■= {A e Cg,e I A(/) > Vp e ^}. 

For an integer t G N anrf G C K[x]2t, taking E = K[x]t, we abbreviate our notation and set 
Ca.t := -CcKMt and Cct.y -CcK^t,^ when K = M. 

4.1. Truncated Hankel operators and radical ideals. In this section, we assume that E 
is a finite dimensional vector space. The following definition for generic elements of Cg,e is 
justified by Theorem 14.61 below. 

Definition 4.4. Let G C {E ■ E) where E is a finite dimensional subspace o/K[x]. An element 
A* G Cg.e is said to be generic if 

(13) lankHfl' = max ranki?f° 

for all subspaces Eq C E. 

If L is a field containing K, we denote by £q ^ :— /^g.b^L, the space obtained by considering 
the vector spaces over L in ((TT|) . We recall here a classical result about generic properties over 
field extensions, which will be used to give a simpler proof of a result that we need from [13] . 

Lemma 4.5. Let K &e a field of characteristic and L a field containing IK. If A* is a generic 
element in ^, then it is generic in £q ^. 

Proof. The space of matrices for A G £q is a vector space spanned by a basis Hi, . . . ,Hi 
over K (resp. L). Let mi, . . . , be new variables and p be the maximal size of a non-zero minor 
G IK[u] of H{u.) :— Y^\^i UiHi. Then for any value of u G IK' (resp. u G L'), the matrix Hixi) 
is of rank < p. Since IK is of characteristic there exists uo G K' with iJ(uo) of rank p, which 
corresponds to a generic element in £q ^ and in £q ^. □ 

Theorem 4.6. Let E be a finite dimensional subspace o/ K[x] and let G d {E ■ E) . Assume 

E 



A* G Cg,e is generic, ie. satisfies US]). Then, Kcr H^, C \/{G) 



A' 

Proof. By Lemma 14.51 A* is a generic element of Cg,e over M or C and thus we can assume 
that IK = K. Let v G V^{G), let 1^ denotes the evaluation at v restricted to {E ■ E) and let 
/ G KeiH^. Our objective is to show that f{v) = 0. Suppose for contradiction that f{v) ^ 0. 

Notice that 1„ and A' := A* + 1„ belong to Lg,e- As A'(/2) = f^{v) 7^ 0, / G Kerfff \ 
KeriJ^, and by the maximality of the rank of iJ^ KeriJ^, ^ KciH^. Hence there exists 
/' G Ker H^, \ Ker Hf. Then, = H^, (/') = i/f (/') + f'{v)l„ implies f'{v) ^ 0. On the other 
hand, 

= Hf, (/')(/) - A' iff ) ^ A{ff) + f{v)f\v) = iJf. (/)(/') + f{v)f'{v) = f{v)nv), 
yielding a contradiction. □ 

4.2. Truncated Hankel operators, positivity and real radical ideals. We first give a 
result which relates the kernel of with the real radical of an ideal (G), when A is positive 
and vanishes on a given set G of polynomials. We start with the following result, which 
motivates our definition of the generic property for a positive linear form. 

Proposition 4.7. For A* G Cg,e,'^, the following assertions are equivalent: 

(i) ranki?^. = maxAg£g ^ ^ rankif^. 

(ii) KeriJf. C Keri/f for 'all A G Cg,e,>- 
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(iii) ranki/^" = maxAeCc e >- rank_ff^° for any subspace Eq C E. 
Call A* G Cg,e,^ generic if it satisfies any of the equivalent conditions (i)-(iii) and set 



^G,E,^ '■= Ker Hj^, for any generic A* e Cg,e,>:- 

Proof (i) (ii): Note that A + A* e Cg^e,^ and Keri/f^^. = KeriJf n Keri?f. (using 
Lemma l3.12p . Hence, rankif^^y^. > rankif^, and thus equality holds. This implies that 
KeriJ|'_|_^, ~ KctH^, is thus contained in KeiH^. 

(ii) =^ (iii): Given Eo C -E, we show that Kerijfp C Kerijf". By Lemma [3?T2| we have 
Kerij|;" C Ker and, by the above, we have Ker C Keri/f . 

The implication (iii) (i) is obvious. □ 
Lemma 4.8. let Go C G C {E ■ E) . Then, /Cgo,b> C /Cg,_e>. 

Proof. Let A S Cg,e,^ be a generic element, so that Ker_ff|' = K,g.e.)^- Obviously, A e 
^Go,E,^, which implies that KerH^ D K.Go.e.)^- O 

Theorem 4.9. Let G C (E-E), where E is a finite dimensional subspace o/R[x]. Then, 
^G,E,)p C \/ (G). 

Proof. Let A be a generic element of Lg.e.)^, so that JCg.e.)^ = Keri7|', and let w e V]r(G'); we 
show that Kei C /(w). As 1,^,, the evaluation at v restricted to (E-E), belongs to £g,e,>^, 
we deduce using Proposition 14.71 that Keri?^ C Ker C I{v). This implies KeriJ^ C 

/(T4(G)) = v^. ~" □ 

Given a subset F C M[x] and t E N, consider for G the prolongation {F \ 2t) of F to degree 
2t, and the subspace E = K[x]f. For simplicity in the notation we set 

(14) ^F,t,:>= := ^(F| 2t>,K[x]t,:>=7 

which is thus contained in \/{F), by Theorem 14.91 The next result (from T^) shows that 
equality holds for t large enough. 

Theorem 4.10. ^ Let F C M[x]. There exists to > such that {JCp^t.-^) = \/{F) for all 
t>to. 

5. Algorithm 

In this section, we describe the new algorithm to compute the (real) radical of an ideal. But 
before, we recall the graded moment matrix approach for computing the real radical developed 
in [13], and the border basis algorithm developed in [T7] . 

5.1. The graded moment matrix algorithm. In the graded approach, the following family 
of spaces is considered: 

= {AeR[x]2t I V/e(F|2t),A(/)=0andVpeM[x]t,A(/)>0}. 

For A e CF,t>. let := i/*'^''. 

Algorithm 15 . II presents the graded moment matrix algorithm described in |12| . 

This algorithm requires in the first step to solve semi-definite programming problems on 
matrices of size the number of all monomials in degree t. This number is growing very quickly 
with the degree when the number of variables is important, which significantly slows down the 
performance of the method when several loops are necessary. The extension to compute the 
radical is also possible with this approach by doubling the variables and by embedding the 
problem over C" in K^". The correctness of the algorithm relies on Theorem 14.101 which comes 
from [ni. 
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Algorithm 5.1: Graded Real Radical 
Input: a finite family F of polynomials of R[x]. 
Set t := 1 and 5 max{deg(/), / G F}] 

(1) Choose a generic A in CF,t>: 

(2) Check wether rankiJ^ = rankff^+^ for some s such that S < s < t: 

(3) If not, increase t := t + 1 and repeat from step (1); 

(4) Compute KeriJ^; 

Output: {/{F) = (kciHl). 



5.2. The border basis algorithm. Algorithm 15.21 presents the border basis algorithm de- 
scribed in [17]. Hereafter, we analyze shortly the different steps. 



Algorithm 5.2: Border Basis 
Input: a family F of polynomials of K[x]. 
Set t:=0,B := {1}, G := and 5 = max{deg(/), / e F}; 

(1) Compute the reduction F of Ft+i on {B)t+i with respect to G; 

(2) Set t' := min{deg(p),p € F,p 7^ 0} - 1 ; 

(3) Compute a minimal G such that (G) := {G+,F) {B+)t'+i; 

(4) Set t" = min{deg(p),p € G n 7^ 0} — 1; Compute B connected to 1 such that 
{B+)t"+i := {B)t"+i ® {G}t"+i; 

(5) Compute a rewriting family G" of Gc+i with respect to Bt"+i] 

(6) If G" 7^ G or S 7^ S or t" < S then set t := t" + 1,B:=B,G := G" and repeat from 
step (1); 

Output: the border basis G of (F) with respect to B. 



In step (1), the reduction of a polynomial p by a rewriting family G for a set B consists 
of the following procedure: For each monomial x" of the support of p which is of the form 
x" = Xi x" x" with x" e S and x" of the smallest possible degree, if there exists an element 
5 = x" — r G G with r e {B), then the monomial x" is replaced by x" r. This is repeated 
until all monomials of the remainder are in B. 

Step (3) consists of the following steps: take the coefficient matrix M = [Afo|Mi] of the 
polynomials in G^ U F where the block A/q is indexed by the monomials in dB'^ and the block 
Ml is indexed by the monomials in B for a given ordering of the monomials, compute a row- 
echelon reduction M of M , and deduce the polynomials of G corresponding to the non-zero 
rows of M. For p € G corresponding to a non-zero row of M, the monomial indexing its 
first non-zero coefficients is denoted 7(p). Notice that (G) := {G'^,F) D {B'^)t'+i contains the 
elements of G+(Gt'). 

Step (4) consists 

• of removing the monomials 7(7?) for p £ {Gt"+i) H {B)t"+i, and 

• of adding the monomials in dB \ {j{p) \ p e G} of degree < t" + 1. 

Step (5) consists of auto-reducing the polynomials p £ G of degree < t" so that j{p) is the 
only term of p in dB. This is done by inverting the coefficient matrix of G with respect to the 
monomials in i9;B. Notice that as {B~^)t"+i ■= {B)t"+i(B {G)t"+i, G is complete in degree t" -|- 1 . 

In step (6), if the test is valid then the loop start again with G a rewriting family of degree 
t with respect to B, which is by definition included in {B^)t- Thus, at each loop, G contains G 
and G+(G) C (G+) n {B+)t+i C (G). 
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The algorithm stops if G" = G and S = S and t > 6. Then t" ^ t, G = G and C+(G) C 
G = G is reduced to by G. If G is a rewriting family complete in degree t for B, we deduce by 
Theorem 1 2 . 1 21 that ttcb is the projection of K[x]t on {B)t along (G 1 1). As B = B, we also have 
t > niax{deg(6) \ b G B} so that G is a border basis with respect to B. As t > 6, the elements 
of F reduce to by G C F. Thus (G) (F). 

It is proved in |17| that this algorithm stops when the ideal (F) is zero-dimensional. Thus 
its output G is the border basis of the ideal (F) with respect to B. 

5.3. K-Radical Border Basis algorithm. Our new radical border basis algorithm can be 
seen as a combination of the graded real radical algorithm and the border basis algorithm. The 
modification of the border basis algorithm consists essentially of generating new elements of 
the (real) radical of the ideal at each loop (step (!') in Algorithm 15. 3p . and to use these new 
relations (which are in the (real) radical by Theorem 14.61 and Theorem 14. 9p in step (3). In the 
case when K = C, a final stage is added to get the generators of the radical of a Gorenstein 
ideal (step (7) below). 



Algorithm 5.3: K-Radical Border Basis 
Input: a family F of polynomials of K[x]. 
Set t = 0, S = {!}, and G = 0; 

(1') Compute a (maximal) S C Bt+i such that S ■ S can be reduced by G onto Bt+i and 
K := GENERICKERNELK(G,i3, 5); 

(1) Compute the reduction F of Ft+i on {B^)t+i with respect to G; 

(2) Set t' := min{deg(|3),p e FUK,p^ O} - 1; 

(3) Compute G such that (G) := {G+,F,K) D {B+)t'+i; 

(4) Compute B connected to 1 and t" < t' maximal such that 
{B+)t-+i := {B)t-+i ® {G)t"+i; 

(5) Compute a rewriting family G" of Gt"+i with respect to B; 

(6) If G" ^GoxB^B or t" < 6 then set t := t" + 1,B:=B,G := G" and repeat from 
step (1); 

(7) i/K = C then [G,B] := Socle(G, 6, A); 

Output: The border basis G of the ideal \/{F) with respect to B. 



The two new ingredients that we describe below are the function GenericKernel (see 
Algorithm [nH]) used to generate new polynomials in the (real) radical, and the function Socle 
(see Algorithm 15. 5p which computes the generators of the radical from the border basis of a 
Gorenstein ideal when IK = C. 

Definition 5.1. Given a rewriting family F with respect to B and S = {x^S • ■ • : x*^' }) 
define F^'°'^ as the following family of polynomials : For all x^' , x^j £ S such that TTp is{'x^^~^^' ) 
exists and is in {S ■ S) , we define k^.+^^(x) ~ x'^'+'^J — 7ri?.B(x'^'+'^J' ) and K^.j^-p. = otherwise. 

With as in Definition 15. li we are going to analyze the corresponding spaces £j?rcd 5, 
Cprad 5 ^, JCprad 5, /C^icd 5 y. Noticc that by construction C {F\t) where t = 2 max{deg(s) | 
s e S} ~ 

The construction of the generic kernel ICprad g (resp., ICprad g^^) is implemented by Algorithm 
15.41 This routine is the one that is executed for finding effectively new equations in the (real) 
radical. 

Notice that primal-dual interior point solver implementing a self dual embedding do return 
such a solution automatically. For a remark on how to use other solvers, see [T^ Remark 4.15]. 
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Algorithm 5.4: GenericKernelk(F, 5) 
Input: A rewriting family F with respect to B allowing reduction for all the monomials 
in S* • 5. 

(1) If K = C, we construct an element A S Cprad g such that has maximal rank, by 
taking a generic clement of the linear space Cpicd g. 

(2) If K = K, we construct an element of A e Cpr^d g^^ such that has maximal rank, 
by computing an element in the relative interior of the feasible region of the following 
semi-definite programming problem: 

- H = {ha,i3)a,pes ^ 

- H satisfies the Hankel constraints Hq^ = 1, ha^js = ha' ^p' if a + /? = a' + /?'. 

- H satisfies the linear constraints haHp^a — for all f3 G S ■ S such that 

(3) Then we compute _fir as a basis of the kernel of H^. 



Output: A family K of polynomials in 




Algorithm 5.5: Socle(G', S, A) 
Input: A border basis G for B connected to 1 and As {B ■ B)* such that is invertible. 

(1) Compute a dual basis of B = {&i, . . . , br} as follows: [di, . . . ,dr] = H^^[bi, . . . , &r] 
where H ~ {A{bibj))i<:ij<r is the matrix of H^; 

(2) Compute A ~ and the matrix iJA = (-'^(^ bj))i<i.j<r by reduction of the 
elements A bi bj by G to linear combinations of elements in B; 

(3) Compute G" = kcr Ha and apply the normal form algorithm to G" U G in order to 
deduce a basis B C B connected to 1 and a border basis G" for B such that 

(G") = (G' U G) = VF. 

Output: A basis B connected to 1 and a border basis G" of for B. 



6. Correctness of the algorithms 

In this section, we analyse separately the correctness of the algorithm over M and C. 

6.1. Correctness for real radical computation. We prove first the correctness of Algorithm 
Oovcr M. 

Lemma 6.1. If G is a rewriting family complete in degree 2t for B, then for A G (G | 2t)^ 

Ker iJ^'^'l' = Ker ijf ' mod (G 1 1) . 
Proof For A e (G | 2 1)-^, we have 

p e Ker Hf""^' ^ A(p g) = Vq G K[x]t 

^ A{bq) = Vg e K[x]t where b £ {B)t = p mod (G 1 1) 
^ A{bb') =OVbe {B)t 
^ 6eKerijf'. 

Therefore 

Ker Hf""^' = Ker iff* mod (G 1 1) , 
which proves the equality of the two kernels modulo (G 1 1). □ 

Lemma 6.2. If G is a rewriting family complete in degree 2t for B, such that K[x]2t = 
{B)2t®{G\2t), then 

{lCG,t> 1 1) = {ICar^d^B^.y 1 1) mod {G\t). 
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Proof. Let A e (G| 2t)-^ be a generic element such that /Cg,*,^ = KerT?^^''''. By Lemma [6.11 
and Proposition I4.7[ we have 

K.G,t,y = KeriJ^'"'' = Kerijf* mod (G 1 1) D JCb,,^ mod (G 1 1). 

Conversely, let A e {G""^)^ be a generic element such that K,G.Bt,h Kcri/f*. As (C'''^) C 
(G I 2t), there exists A G {G\2t)^ which extends A to K[x]t. Then we have 

/CcB^t = KerHf' =Kcriff 

= Keri/J'""!* mod (G 1 1) D /Cg,*,>: mod (G | i). 

□ 

Lemma 6.3. If Algorithm ] 5.^ terminates with outputs G andB, then (G) = \/{F) and B is a 
basis o/K[x]/y(F). 

Proof. If the algorithm stops, all boundary polynomials of G+(G) reduce to by G. By Theorem 
[2T21 for all t we have K[x]2t (6)2* ® {G\2t). As /Ccrcde^^ = {0} by Lemmai^l we deduce 
that 

JCG,t> C (G|t). 
By Theorem 14. 10[ there exists sq such that 

where / = (F). By lemma for t > sq, 

ICF,so.y c ICG,t,y c (G|t) c v^, 

which implies that (G) = v^. □ 

Proposition 6.4. Assume that Vr(-F) is finite. Then the alaorithm \5.3\ terminates. It outputs 
a border basis G for B connected to 1, such that K[x] = (B) (B (G) and (G) = \fl . 

Proof. First, wc are going to prove by contradiction that when the number of real roots is finite, 
the algorithm terminates. 

Suppose that the loop goes for ever. Notice that at each step cither G is extended by adding 
new linearly independent polynomials or it moves to degree t + 1. Since the number of linearly 
independent polynomials added to G in degree < t \s, finite, there is a step in the loop from 
which G is not modified any more. In this case, all boundary G-polynomials of elements of G 
of degree < t are reduced to by Gt. By Theorem l2.12[ we have 

(15) U.[At = {B)t®{Gt\t). 

We have assumed that the loop goes for ever, thus this property is true for any degree t. By 
Theorem 14. 10[ there exists sq such that 

As any element of {F \ sq) reduces to by the rewriting family Ggg , we have {F \ sq) C {Ggg \ sq). 
By Lemma l48l we deduce that 

I^F.,So/2.t C ICG,^,so/2,y- 

For a high enough number of loops, the set Gsq is not modified and wc have ^Ga^.s^^/j.^ " {'^i- 
Applying Lemma 16.21 using Equation (|15l) . we have 

K^Gso,so/2,y C (Gso I So) 

By construction Ggg C v^, thus 
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Let Bo C M[x] which defines a basis in M[x]/ \^ and of smallest possible degree and let do be 
the maximum degree of its elements. Then any monomial m of degree + 1 is equal modulo 
(,+1 to an clement b in (Bq) of degree < do- 
By Theorem Em 

(Gdo+i I c?o + 1) = ( v^)do+i> 
thus m — b £ {Gdo+i | do + 1) so that any monomial of degree do + 1 can be reduced by G to a 
polynomial in K[x]rf„. Thus B = n/gG(7(/))^ C 1R[x](Jq is finite and the algorithm terminates. 

By Lemma [Ol the algorithm outputs a border basis G with respect to B connected to 1, 
such that (G) ^ Vl- □ 

6.2. Correctness for the radical computation. In this section, we show the correctness of 
the algorithm for radical computation, that is with K = C. 

Proposition 6.5. Assume that Vc{F) is finite. Then the algorithm 1 5. 3\ terminates and outputs 
a border basis G for B connected to 1, such that (G) = \/l and C[x] ~ {B) © 

Proof. Since the family G contains the polynomials constructed by the normal form algorithm 
[17j and as Vc {!) is zero-dimensional, the normal form algorithm terminates and so do algorithm 
15.31 When the loop stops, all boundary polynomials of G~^{G) for any degree reduce to 
by G and K.Qiod ^ ~ {0}. By Theorem I2.13[ G is a border basis with respect to B. Let 
A e {B ■ B)* such that /Cgrcd g = KeriJ^. By definition of A and normal form property, if 
/ G {B ■ B) n (G) then / reduces to by G and A(/) = 0. This shows that we can extend A 
to A e C[x]* by A = A on (B) and A = on (G). We deduce that (G) = Keri?^ and that 
A A = C![x]/ Keri/^ = C[x]/(G) is Gorenstein. Let di, . . . , d^ be the dual basis of B for Q\ and 
A ^ X]i=i bidi. By Theorem 13.101 Kei H^,^ computed in the function Socle, yields a new 
basis B' connected to 1 and a new border basis G' such that (G') = a/T- Q 

7. Examples 

This section contains two very simple examples which illustrate the effect of the SDP solution 
in one loop of the Real Radical Border Basis algorithm. The results in the next example are 
coming from a C++ implementation available in the package newmac of the project mathemagix. 
It uses a version of lapack with templated coefficients and sdpa0 with extended precision so 
that all the computation can be run with extended precision arithmetic. 

7.1. A univariate example. Wc give here a simple example in one variable to show how the 
real roots can be separated from the complex roots, using this algorithm. We consider the 
polynomial f = x'^ — x^ — x + 1 with a single real root a; = 1 of multiplicity 2. In the routine 
GenericKernel of the algorithm, a 4 x 4 matrix H is constructed and the linear constraints 
deduced from the relations x* = x^ + x — 1, x^ = x^ + x'^ — 1, x^ = 2x^ — 1 modulo / imposes 
the following form: 

/la b ^ \ 

^ _ a b c c+a— 1 

6 c c + a — lc + 6— 1 
\ c c+a-1 c+b-1 2c-l / 

where a = A{x), b = A(a;^), c = A{x^). The SDP solver yields the solution 



/I 


1 


1 


1\ 


1 


1 


1 


1 


1 


1 


1 


1 


V 1 


1 


1 


1 / 



■http : //sdpa. sourcef orge .net/ 
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which kernel is {x — l,x'^ ~ x,x^ ~ x"^). Thus the output of the algorithm is {x — 1) the real 
radical of (/), the basis B = {1} and the real root x = 1. 



7.2. A very simple bivariate example. Let fi ^ x'^ + and F = {/i} C ]R.[x,?/]. The 
algorithm computes the following: 

— We compute GenericKernel in degree 1 by choosing S = {l,x,y} with S ■ S D 
support /i . 

The SDP problem to solve reads as follows: find h = [a, b, c, d] £ M"* such that 



H 



1 a b 

a c d 
b d -c 



^0 



and has of maximal rank. Here a = A{x),b = A{y),c ~ A(x^),(i — A{xy). The 
condition H )p Q implies that 

- c = 0, 

- a = 0,& = 0,d = 0. 

and consequently that keri/ = {x, y). Thus x, y are returned by GenericKernel and 
added to F. 

— After one iteration the border basis algorithm stops and we obtain B = {1} and 
V(x2+J/2) = (a;,y). 



7.3. Numerical example. The tables below compare the size of the SDP problems to solve in 
our approach and in the method described in |12j . The degree indicates the degree in the loop 
of the Border Basis Real radical algorithm, n.sdp is the size of matrices in the corresponding 
SDP problem and n. constraints the number of linear constraints involved, t is the degree of 
the relaxation problem in |12) and n.sdp grad. rel. the size of matrices in the corresponding 
SDP problem. 



Katsura 4 


degree 


n.sdp 


n.constraints 


t 


n.sdp grad. rel. 


2 


5 


5 


2 


56 


4 


11 


67 


2 


56 


6 


16 


176 


2 


56 



Katsura 5 


degree 


n.sdp 


n.constraints 


t 


n.sdp grad. rel. 


2 


6 


6 


3 


84 


4 


16 


146 


3 


84 


6 


26 


479 


3 


84 
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bifur 


degree 


n.sdp 


n.constraints 


t 


n.sdp grad. rel. 


2 


4 


2 


8 


165 


4 


9 


32 


8 


165 


6 


16 


150 


8 


165 


8 


25 


446 


8 


165 


8 


16 


152 


8 


165 


8 


16 


153 


8 


165 


6 


16 


158 


8 


165 





16 


lOZ 


o 

8 


165 


4 


9 


34 


8 


165 


6 


16 


168 


8 


165 


6 


16 


169 


8 


165 


4 


9 


36 


8 


165 


6 


16 


177 


8 


165 


4 


4 


3 


8 


165 


4 


8 


37 


8 


165 



The tables below give the time for computing the real radical with the solvers sdp4j or csd] 
integrated into the border basis algorithm available in the package newmac of mathemagix. 



Example 


Tk 


Gen. Ker. 


CSDP 


SVD Drop 


Deg 


Degc 


Nr 


Nc 


Tc 


Precision 90 


katA 


22.479s 


22.281s 


22.1645 


le - 12 


4 


4 


12 


16 


0.06s 


katb 


146.49s 


146.29s 


145.64s 


le - 12 


5 


5 


16 


32 


0.165s 


cyclo 


10.839s 


10.7646s 


10.6243s 


le-20 


5 


5 


4 


16 


0.03s 


robot 


41.84s 


41.52s 


41.26s 


le - 19 


6 


8 


4 


40 


1.3s 


Precision 120 


katA 


22.557s 


22.28s 


22.16 


le - 14 


4 


4 


12 


16 


0.06s 


katb 


146.59s 


146.39s 


145.1s 


le - 12 


5 


5 


16 


32 


0.17s 


cyclo 


10.839s 


10.7646s 


10.6243s 


le-20 


5 


5 


4 


16 


0.03s 


robot 


42.884s 


42.5216s 


42.2447s 


le - 19 


6 


8 


4 


40 


1.4s 



A precision of 90 or 120 bits is used during the computation but unfortunately the SDP solver 
is very, very, very slow for this precision. A strange behavior/bug of the parameter used in the 
relaxation of the barrier function is observed. The solution of this problem is in progress. 



Example 


Tk 


Gen. Ker. 


SDP A - GMP 


SVD Drop 


Deg 


Degc 


Nr 


A^c 


Tc 


Precision 90 


katA 


4.18s 


4.12s 


3.26 


le - 18 


A 


A 


12 


16 


0.06 


katb 


26.28s 


26.01s 


23.16s 


le - 18 


5 


5 


16 


32 


0.165 


cyclo 


10,95 


10.77 


10.64 


le-20 


5 


5 


4 


16 


0.03 


robot 


19,84 


19.52 


19.26 


le - 19 


6 


8 


4 


40 


1.3s 



Using SDPA-gmp as the solver allows us a great improvement in efficiency though we expect 
futhcr improvements improving both the way connection with SDPA-gmp is operated and better 
tuning the parameters SDPA. 
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